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Abstract. We review and study the roles of quantum and classical fluctuations 
in recent cavity-optomechanical experiments which have now reached the quantum 
regime (mechanical phonon occupancy < 1) using resolved sideband laser cooling. 
In particular, both the laser noise heating of the mechanical resonator and the form 
of the optically transduced mechanical spectra, modified by quantum and classical 
laser noise squashing, are derived under various measurement conditions. Using this 
theory we analyze recent ground-state laser cooling and motional sideband asymmetry 
experiments with nanoscale optomechanical crystal resonators. 
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1. Introduction 

Mechanical resonators provide an interesting and useful system for the study of quantum 
theory in the mesoscopic and macroscopic scales [H El El IH El |6]. Technological 
progress in nano- and microfabrication have made coupling of motion to a wide variety 
of systems possible, and experiments demonstrating mechanical resonators coupled 
to optical cavities [fl El ® Hffl HD H21 H31, microwave resonators [H] 115] HE] 117]. 
superconducting qubits [HI [191 [20] , c °ld atoms [21], and defect centers in diamond [22] 
have been forthcoming. These advances have led to the cooling of mechanical systems 
to their quantum ground states [121 dHl [23], and the observation of nonclassical 
behaviour [191 1231 [25] . Of all the systems to which mechanical systems have been shown 
to couple, it has long been realized [JJ that light has many distinct advantages. Most 
importantly light sources can be made quantum limited (lasers) and lack thermal noise 
even at room temperature. This allows, for example, to cool optomechanical systems 
close to their ground states using radiation pressure forces, and recently experiments 
starting with modest cryogenic pre-cooling [HI [231 EH EH], have been successful at 
achieving this task. 

Considering the importance of cooling and thermometry in cavity-optomechanical 
experiments, it is important to understand the noise processes which can lead to both 
heating and systematic errors in thermometry. In this paper we focus on analyzing these 
effects theoretically, and compare to recently performed experiments with nanoscale 
optomechanical crystal devices. In Section [2] a theoretical treatment of optomechanical 
cooling and thermometry particularly suited for understanding the propagation of 
noise is given. The regime of operation analyzed is that relevant to our recent 
experiments [231 123]. i.e. the driven weak-coupling, sideband-resolved regime where the 
cavity decay rate is larger than all other rates in the system except for the mechanical 
frequency. In this regime, the standard input-ouput formalism [271 EH EH] for analyzing 
the linearized system is applied. Two different methods of thermometry used in recent 



experiments [23J E3J are treated in sections 2.2 and 2.3, respectively. The former, 
involves directly measuring the light scattered by mechanical motion and calibrating 
its intensity given a set of system parameters, while the latter requires comparing the 
emission and absorption rates of phonons from mechanical subsystem, and observing 
a non-classical asymmetry analogous to that seen in much earlier experiments with 
trapped atoms and ions [301 EU [32]. Both methods are susceptible to laser phase noise, 



which is discussed in section 2.4, through both noise induced heating [33], and systematic 
errors in thermometry caused by noise squashing [33] US] and anti- squashing. Finally, 
in Section [3] we review recent laser cooling and sideband asymmetry measurements 
of optomechanical crystal cavities near the quantum ground state of their mechanical 



motion, and compare these results with the measured phase noise (see Appendix E) of 
the external cavity semiconductor diode lasers used in these experiements. 
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2. Theory 
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Figure 1. a, Bi-directional evanescent coupling geometry, in which the "transmitted" 
field goes into the forward a ou t waveguide channel and the "reflected" field goes into the 
backward a' out waveguide channel. This is the coupling geometry we will be focusing 
on in this work, in which the laser input channel is aj n and the detection channel is the 
forward waveguide channel, a ou t- b, Double-sided end-fire coupling geometry. This is 
the geometry one would have in a Fabry-Perot cavity. Note that in this geometry what 
we call the "transmission" and "reflection" channels are typically opposite of that in 
the evanescent coupling geometry (a direct map between the two geometries would 
relate the "reflected" channel in (a) to the conventional transmission of a Fabry-Perot, 
for instance), c, Single-sided end-fire couplng geometry. This is the ideal measurement 
geometry, in which all of the optical signal that is coupled into the cavity can be, in 
principle, collected and detected in the a out channel. In principle one does not have 
excess vacuum noise coupled into the optical cavity, only that from the laser input 
channel a m . 



The optomechanical interaction between a mechanical system and an optical field 
occurs through radiation pressure, a force proportional to the optical field intensity. 
This can be modeled by a Hamiltonian 

H = huj a)a + hu m o&b + HgaJa(b> + b), (1) 

with a and b the annihilation operators for photons and phonons in the system. In 
the presence of a laser emitting light at frequency ul, it is convenient to work in 
an interaction frame where u Q — > A in the above Hamiltonian with A = uj Q — ul- 
To incorporate the effect of the environment, we use the quantum-optical Langevin 
equations for the system [27J E3 129] . 

b(t) = - (iu m0 + & - igrfa ~ 

and 

a(t) = - (iA + |) a - iga(& + b) - ^J2a hi (t) - V^a^t). (2) 

Here ai n (t) and a- m ^(t) are the quantum noise operators associated with extrinsic 
(input /output) and intrinsic (undetected) optical loss channels, respectively. Here we 
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assume bi-directional evanescent waveguide coupling (see Fig. [T]) to the optical cavity in 
which the total extrinsic cavity loss rate into both directions of the coupling waveguide 
channel is K e , with the uni-directional input coupling rate being half that at K e /2. The 
total optical cavity (energy) decay rate is given by k, with Kj — Kj Kip /2 denoting all the 
optical loss channels which go undetected. [|] The noise operator b- m (t) arises from the 
coupling of the mechanical system to the surrounding bath degrees of freedom, which 
in most current systems resides in a high temperature thermal state with average bath 
occupancy n& » 1. 

We linearize the equations about a large optical field intensity by displacing 
a — > «o + a - Such an approximation is valid for systems such as ours where go « k, 
and the optical vacuum alone only marginally affects the dynamics of the system, i.e. 
the vacuum weak coupling regime. All experimental systems to date are in this regime. 
For systems in the vacuum strong coupling (go > k) regime, more elaborate treatments 
taking into account the quantum nature of the nonlinearity must be pursued [35j EH1 EZ] . 
In the Fourier domain the operators for the mechanical and optical modes are found to 
be 

-■sffib in (u) iG(a(uj) + a 1 " (a;)) 



and 



b(u) 



aiu) 



i(u m0 -co) + 7i/2 



K e /2a,- u 



i(u m0 7i/2 

iG(b(u) + ft(u)) 



k a 



m.i 



i(A + k/2 
respectively, where G = g\a \. 

Using equations ([3]|4]) we arrive at the operator for the mechanical fluctuations 



(3) 



(4) 



b(u) 



+ 



l[OJ m 

iG 



+ - 



- u) + 7/2 

V ' K e /2a- m {uj) + ^K'di^u) 

i(A — to) + k/2 i(oo m — w) + 7/2 

iG V^MnM + V^aL(w) 



(5) 



-i(A + lu) + k/2 i(u m -u) + 7/2 
where u m = u m o + 5u m is the optical spring shifted mechanical frequency and 7 = 
7i + 7om is the optically damped (or amplified) mechanical loss-rate. Expressions for 
the optical springing and damping terms are given by 

1 1 



and 



7om 



|GfIm 



2|G| 2 Re 



z(A 



to. 



,) + k/2 -i(A + Lu m ) + k/2 



z(A 



t ) + k/2 -i(A + w m ) + k/2 



(6) 



(7) 



I In an idealized measurement, all photons lost by the optical cavity are lost to the detected channel 
so that k' = 0. In this ideal case, satisfied by the single-sided end-fire coupling geometry of Fig. [lj;, 
the vacuum fluctuations from ports other than the detector port never enter the optical cavity. For 
a bi-directional coupling scheme we have k' > K e /2, and there is information lost in the backwards 
waveguide direction about the mechanical state. 
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respectively. From these expressions it is evident that in the sideband-resolved regime 
the maximum optical damping occurs for a laser red-detuned from the optical cavity 
with A = u m , resulting in a damping rate of 7om = 4|G| 2 /k. The ratio between this 
optical contribution to the mechanical damping and the intrinsic mechanical damping 
is the co-operativity, C = 7om/7i- 



2.1. Quantum-limited laser cooling and damping 

The expression for the noise power spectrum of the laser driven mechanical system 
can be calculated using eqn. ^ for b(u). More specifically we calculate Sbb(u) (see 
Appendix B), corresponding in the high-Q regime to the ability of the mechanical 
system to emit noise power into its environment |B|. The area under S^u) is the 
average mode occupancy of the mechanical quantum oscillator. In the absence of the 



optical coupling to the mechanics (G = 0), the result in Appendix B is obtained. 
Allowing for optical coupling and including the optical noise terms, we arrive at an 
expression involving the correlations (a-" n (u;)di n (u/)), (aj ni (u;)a irii i(a/)), {di Q (u)d\ n (uj')) 
and (a inti (u)a\ ni (u)')) that must be calculated from the properties of the optical bath. 
Assuming that our source of light is a pure coherent tone, and thus the optical bath is 
in a vacuum state, as is approximately the case in many optical experiments, the former 
two correlations can be set to zero, while the latter two give 5(u> + u>'). As described 
in Appendix B for the mechanical system which is in contact with a thermal bath of 
occupancy n&, we have noise input correlations of {b in (u)bl n (u')) = (n& + l)6(u + u/) and 
(bl n (oo)b- m (oo')) = n b 5(oo + u/). The expression for Sbb{u) is then found to be, 

(u m + u) 2 + (7/2) 2 

where rif(u), the back-action modified phonon occupation number, is given by 

liUb \G\ 2 k 1 

7 7 (A - u) 2 + 

In the driven weak- coupling regime (re » 7), the mechanical lineshape is not strongly 
modified from that of a Lorentzian, and nj(u) can simply be replaced by nf(—u m ) 
in eqn. Q. An input laser beam tuned a mechanical frequency red of the cavity for 
optimal laser cooling (A = w m ), results in a back-action modified average mechanical 
mode occupation number equal to 

/~\l linb . ^ OM ( K V fin\ 



n f (u) = — + —^—^ — 2 , rzw - ( 9 ) 



7 7 \ 

The term n q bi = (K/4u m ) 2 is the quantum limit on back-action cooling, as derived in 
[38| 139] using master-equation methods and in [ID] by taking into account the spectral 
density of the optical back-action force. This small (in the good cavity limit) residual 
heating comes from the non-resonant scattering of red pump photons, to one mechanical 
frequency lower, or a total of 2u m detuned from the optical cavity. 
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We note briefly that for other laser detunings different back-action occupancies are 
achieved, such as 

4IGI 2 / k x2 



and 



<n>| A=0 = n 6 + ^ — (11) 



<n> A=-w m = + , (12) 

7 7 

where again the resolved-sideband limit is assumed, and for A = —u m one has 
amplification of the mechanical motion with 7qm = ~ 4:\G\ 2 /k. 

2.2. Thermometry with the cooling beam 

One of the simplest methods of inferring the mechanical mode occupancy is to detect 
the imprinted mechanical motion on the cooling laser beam itself. Upon transmission 
through the cavity-optomechanical system, the laser cooling beam, typically detuned to 
A = uj m in the resolved sideband regime, preferentially picks up a blue-shifted sideband 
at frequency u>l + w m (^ w B ) due to removal of phonon quanta from the mechanical 
resonator (anti-Stokes scattering). Upon detection with a photodetector, the beating of 
the anti-Stokes sideband with the intense cooling tone produces an electrical signal at 
the mechanical frequency [TOj SH [23] . By careful calibration and accurate measurement 
of the magnitude of this signal, and through independent measurements of other system 
parameters such as g, \cxq\, K e , and k, the mechanical resonator's average phonon number 
occupancy can be inferred. 

The optical fluctuations in the transmitted laser cooling beam at the output port 
of the optomechanical cavity are given approximately in the sideband resolved regime 
by 

OoutMiA=u) m » t(u; A)a in (u) + n opt (uj; A)a in>i (u;) + s 12 (uj; A)b in (u), (13) 
where t, n opt , and Si 2 are the scattering matrix elements evaluated for a laser cooling 



beam of red-detuning A = uj m (see Appendix D). This expression is derived using eqns 



(|3j4) and input-output boundary condition a out (u) = a in (uj) + K e /2a(u). Expressions 
of this form have been used previously to analyze the propagation of light and sound 
through an optomechanical cavity in the context of state transfer [121 H3] . In this case 



we have simplified eqn. (13) by ignoring input noise terms from the creation operators 
a\ n (u) and a\ ni (u). These terms gives rise to the quantum-limit of laser cooling found 
in eqn. (10) above, but insignificantly modify the optically transduced signal of the 
mechanical motion as long as (n) » n q u- 

The strong cooling laser tone beats with the optical noise sidebands, generating a 
photocurrent proportional to a out (t) + a|, ut (t), 

I{u)\A=w m = t(u)d in (u) + n op t(w)a in> i(w) + s 12 (u)b iQ (u) + 
+ h.c.(-u) 
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where h.c.(— u) is a convenient short-hand (f(oj) + (f(—u)y = f{oo) + h.c.(— u)). The 
resulting photocurrent power spectral density as read out from a spectrum analyzer is 
given by 



S n (u) 



rCO 

du' (P(u)i(u')) 

«Mf(- u r(«ioM«L(w')) 



+ n opt (u)n opt (-u'y(a inii (u)a\ n -(uj')) 
+ s 12 (u)s 21 (-u'y(b in (u})bl(u')y 

+ s 12 (-u) S21 (ur(HMhn(u')>du', (14) 

where we have assumed the same optical (vacuum) and mechanical (thermal) noise 
correlations as above in evaluating eqn. fl8J). Using the normalization property of the 
scattering matrix coefficients (|t(u;)| 2 + |^o P t(w)| 2 + |si 2 (u;) | 2 = 1), we find the simplified 
expression 

Sn(co)\ A ^ m = 1 + n b (\s 12 (u)\ 2 + M~^)| 2 ). (15) 
Substituting for the expression of the phonon-photon scattering element s\ 2 (u) given in 



Appendix D for A = u m yields 



Su(u)\A=„ m = 1 + ^^£ 66 (o;;<n», (16) 
2k k 

for the transduced noise power spectral density, where (n) is the actual mechanical 
mode occupancy including back-action effects of the cooling laser (see eqn. (10)). 

Several points are worth mentioning regarding this expression. Firstly, the signal 
to noise goes as the coupling efficiency r] = k c /2k. Secondly, the detected signal is 
proportional to (b^b) as opposed to (b^b) + 1/2, and so the resulting signal is exactly 
what would be expected classically, vanishing as the temperature and phonon occupation 
go to zero. In other words, this measurement is insensitive to the zero-point motion of 
the resonator. The spectral density Sbb(u,(n}), represents the ability of the mechanical 
system to emit energy [B]. By tuning the laser to A = u m in the sideband- resolved 
regime it is exceedingly unlikely for the tone to drive the mechanics (through Stokes 
scattering), and so we gain little information about how the mechanical system absorbs 
energy from the optical bath. Finally, we note that equation (16) is general and holds 
for both low and high cooperativity. 



2.2.1. Intepretation as quantum noise squashing Though the scattering matrix 
formulation provides a consistent and systematic way of deriving the form of the detected 
signals, it does so by elimination of the position operator from the equations. It 
is interesting to reinterpret the experiment as a measurement of the position of the 
mechanical system [UJ, and we attempt to do so here^J For simplicity, the perfect 

§ A much more thorough treatment of the implications of quantum back-action and measurement 
theory in this type of system is presented in a recent work by Khalili, et al. [33] 
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coupling condition is assumed, i.e. K e /2 = k. The output signal is then given by 



The normalized heterodyne current is found to be 

2G 

I(t)\A=. m = -ia- m (t) + ial Q (t) + —(b(t) + 6t(t)) (17) 



and so it would seem that the signal is composed of optical shot-noise and a 
component which is proportional to x, making Sjj = 1 + const x <x 2 ). This however 
contradicts the above derivations which show that Su = 1 + const x (n) for A = u m . The 
inconsistency comes after careful calculation of the correlation function (/(£ + r)J(i)). 
In fact, a\ n (t) and b(t) are correlated, and the view that the shot-noise simply creates a 



constant noise floor is incorrect. Proper accounting for the correlations (see Appendix 



Cj) leads us again to eqn. (16), showing that the measured quantity is Sbb, and the area 
of the detected spectrum is proportional to (n). 

The blue-side driving with A = — uj m causes the opposite effect, i.e. quantum noise 
anti-squashing. The squashing and anti-squashing are signatures of quantum back- 
action. This effect is in spirit similar to classical noise squashing which we study in 



Section 2.4.2 , where correlations between the noise-induced motion and classical noise of 
the detection beam destructively interfere at the photodetector. It is important to note 
that this signature of quantum back-action does not involve detection of quantum back- 
action heating, and can be apparent at arbitrarly low powers, far below that required 
to reach the standard quantum limit. 



2.3. Motional sideband asymmetry thermometry 

An alternate method of measuring the temperature of the mechanical subsystem, one 
which uses the mechanical zero-point motion to self-calibrate the measured phonon 
occupancy, involves comparing the measured signal from a weak probe beam (low 
cooperativity) at both A = +u m in the sideband resolved regime [21]. In such 
experiments, the mechanics can be either laser cooled with a different laser and/or 
optical cavity mode, or the system can be cryogenically pre-cooled to a temperature 
which requires no further cooling to approach the quantum ground state. As the 
optical read-out beam can be arbitrarily weak in such measurements, it only marginally 
affects the dynamics of the mechanical system [15] [11] [21]. By working at low read- 
out beam power, such that the optically-induced damping and amplification rates are 
much smaller than the bare mechanical linewidth, optical back-action by the probe 
beam only minimally affects the dynamics of the mechanical system and measurements 
can be taken at detunings both red (A = u m ) and blue (A = — u m ) of the cavity 
without triggering any optomechanical instabilities |3B]. Operating in the resolved 
sideband regime allows for the separate cavity filtering of the Stokes and anti-Stokes 
motionally induced sidebands on the probe beam, which are respectively proportional 
to (n)+l and (n). It can be shown that the additional vacuum contribution to the Stokes 
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scattering, which provides the intrinsic calibration for (n), arises in these measurements 
equally from the shot noise on the probe laser and zero-point motion of the mechanical 
resonator [H]. We will also see in the following sections, that such a measurement at 
both A = +u m can provide additional resilience to systematic errors from non-idealities 
such as laser phase noise. 

We derive here the blue-detuned (A = — uj m ) result analogous to the red-detuned 
(A = ui m ) laser cooling case given above in eqn.(Jl6|). In the sideband- resolved regime, 



the approximations that led to eqn. (13), lead to a similar expression in the case of 
A = — bj m for the electromagnetic field output from the optomechanical cavity 

aout(w)| A =- Wm » t(u; A)a in (u) + n opt (u; A)a- m;i (u) + s 12 (u; A)foJ n (u;). (18) 

where we have neglected the terms proportional to the photon noise creation operators as 
their effect is again minimal on the optically transduced signal of the mechanical motion. 
Such a scattering relation, whose exact form is shown in the Appendix D allows the 
optomechanical system to act as an amplifier, and has been studied experimentally at 
microwave [17] and optical frequencies [37], and studied more generally in the context 
of optomechanics by Botter, et al. [H]. The scattering elements satisfy the equation 
|i(o;)| 2 + |ra opt ((jj) | 2 + |si2(w)| 2 = 1, which along with the standard bath correlation 
relations used above, allows us to write 

S , //(w)| A =- Wm = 1 + K + l)(|ai 2 (w; A)| 2 + |s 12 (-w; A)| 2 ) (19) 

= l + ^^Wu;;<n», (20) 

where (n) is the actual mode occupancy including back-action of the laser input (see 
eqn. (Jl2|)). As before, the signal lies on top of a flat shot noise background of unity, and is 
proportional to the detection efficiency rj and the measurement rate 7om- Now, however, 
the signal is proportional to the creation operator spectral density S b i b t(u;(h)), which 
itself is proportional to (n) + 1. The spectral density S b i b i(u;(h)) can be interpreted 
as the mechanical system's ability to absorb energy [5], which even at zero temperature 
(occupation) can absorb energy through spontaneous scattering process which arises 
due to the zero-point motion of the mechanical resonator. 

For a constant laser driving power the optomechanical damping and amplification 
rates for detuning A = +co m are equal in magnitude but opposite in sign, with 
7+ = 7j + |7om|, where |7om| = 4|G| 2 /k in the sideband resolved, weak coupling 
regime. Weak probing entails using a probe intensity such that |7om| « 7i, or C r « 1, 
where we define C r = |7om|/7i as the read-out beam cooperativity. In this limit the 
mechanical mode occupation numbers for A = +u m detunings are given approximately 
by (n) + = 7jnb/7+, where n b is the mechanical mode occupancy in absence of the 
probe fielc||] Denoting the integrated area under the Lorentzians in eqns. (16) and 
(20) as I + and respectively, we find a relation between their ratios and the read- 
out cooperativity which provides a quantum calibration of the unperturbed thermal 



Referring to eqns. (10) and (12), this is an accurate relation if C r « rib- 
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occupancy [2 



_ J -/ J + 11 
V = 1 + C r ~ 1 - C r ~ r7 b (21) 

2.4- Laser phase noise 

Although various other noise sources (laser intensity noise, internal cavity noise [15], etc.) 
can be treated similarly, here we focus on laser phase noise as it is typically the most 
important source of nonideality in the laser cooling and thermometry of cavity opto- 
mechanical systems. The effect of phase noise on optomechanical systems has already 
been studied at great depth in the context of heating [33] of the mechanical resonator 
and entanglement [1H1 EDI EI] of light and mechanics. However, laser light often acts 
as both a means by which the mechanical system is cooled as well as its temperature 
measured, and thus laser noise can effect both the true and inferred mechanical mode 
occupancy. Here we complement previous studies of laser noise heating with a unified 
analysis that also quantifies the effects of quantum and classical (phase) laser noise on 
the optically-transduced mechanical mode spectra. 

The optical laser field amplitude input to the optomechanical system, in a rotating 
reference frame at frequency oo E and in units of -y/photons/s, we denote by E Q . Due 
to processes internal to the laser, some fundamental in nature, others technical, this 
amplitude undergoes random phase fluctuations which are captured by adding a random 
rotating phase factor [52] 

E (t) = \Eo\e«®. (22) 

As long as the phase fluctuations are small, we expand this expression to first order 
yielding E (t) w \E \(1 + Uf>(t) + 0((j) 2 )) [33]. Then 

<^(r) J E o (O)>=|E o r(l + <0(r)0(O)». 

In this way, we can express the noise power spectral density of the optical field amplitude, 

S EE (u), as 

S EE (oo) = \E \ 2 (2ti5(uj) + S H (u)), (23) 

where we've also used the realness of <j>{t) to set S^oS) = S^{u). 

This relates the phase noise power spectral density to the optical power spectum 
of the noisy laser beam, with the optical power away from the carrier at u = due to 
phase noise. This phase noise can then be taken into account as an additional noise 
input to the cavity, 

ai n ,totM = «inM + flin^M, (24) 

where a-^^u) is a stochastic input with (a\ n ^(a;)ai n ^(c</)) = S EE (u)6(u + u/) (here the 
averages used for correlation functions correspond to classical ensemble averages and 

There is, however, an additional subtlety when performing mechanical mode 
thermometry with a laser beam affected by phase noise; correlations between the positive 
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and negative frequency components of the phase noise can cause cancellations in the 
optically transduced signal, and therefore must be carefully taken into account. For 
example, for a pure sinusoidal tone phase modulated onto a laser we have, 
Eqc^^ ~ Eq (e^ c COSUJt+i @ s sinw *^ 

E (l + 1 -{(3 S + i(3 c )e wt 1 -{(3 S - ifoe** 

The positive and negative frequency optical sideband amplitudes are negative complex 
conjugates of one another. More generally, the positive and negative frequency 
components of the noisy optical input have the following relation for an optical signal 
with phase noise, 

= ( 25 ) 

where ct[~^(w) = 6(—uj)a in ^(uj), and a[^(uj) = 6(u)ai n ^(u). The total phase noise signal 
can then expressed as a- m ^(u) = aj^(u>) — a\^(—u>). For calculations that follow, this 
explicit separation of positive and negative frequency phase noise components is useful 
in simplifying calculations of the optically transduced mechanical motion. In terms of 
positive frequency phase noise components only then, we have 

<4iiV)O w/ )> = S EE (u)5(u + u')9(u), (26) 

with a similar relation holding for the negative frequency components of the phase noise 
input and the negative frequency optical power spectrum. 

2.4-1- Heating To find the actual thermal occupation of the mechanical system in the 
presence of phase noise on the laser cooling beam we use once again eqn. ([5| for b(u), 
replacing a- m (cu) with d inttot (uj) which includes the classical phase noise on the input laser. 
From the non-zero correlation (a\ tot (u>)di nt tot(^>')} = S E e(u)5(uj + u') for > 0, we 
find another source of noise phonons, in addition to those coming from the thermal bath 
and quantum back-action of the laser light. This is expressed as new terms proportional 
to S EE (u) in the noise spectrum of the mechanical motion given by eqn. 



where 



[u m + uj) 2 + (7/2) 2 

jiU b \G\ 2 k 1 + (K e /2n)S EE (uj) 
nf,<p{u)\ A = Urri = — - + 



7 7 (A-lu) 2 + (k/2) 2 
\G\ 2 k {kJ2k)S ee {u) 
7 (A + lu) 2 + {k/2) 2 ' 1 ' 

As before, assuming a high mechanical Q-factor and the driven weak-coupling regime 
(k » 7), we substitute nf^{— ui m ) for rif^uS) and relate it to the average mode occupancy 
in the presence of laser phase noise, (n)^, 

2 



1 li^b , 7om 



rve 



4oo m J \2k 



(29) 
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where we have defined = SEE^m) [33J. 

The additional phase noise heating in equation (29) can be understood as the 
product of the number of noise photons present in the light field at a mechanical 
frequency detuned from the central laser frequency (n$), multiplied by the efficiency with 
which they are coupled into the cavity (k £ /2k), and finally multiplied by the efficiency 
with which the light field couples to the mechanics (7om/t)- For detuning A = u m used 
in resolved sideband laser cooling, the optomechanical system only samples the input 
laser phase noise at a mechanical frequency blue of the central laser frequency, and thus 
the relationship between the negative and positive frequency components of the laser 
phase noise has no role to play in heating in the sideband-resolved regime (as we will 
see below, this is not the case in the thermometry). 




Figure 2. a, and b show the spectral density of the detected signal, Sh(uj) for a large 
phonon occupation n » 1, about the mechanical frequency. The contributions from 
phase noise (in orange), thermal brownian motion (in blue) and shot-noise (in grey), 
are highlighted. In a, the effect of noise squashing, arising from the cancellation of 
phase noise near the mechanical frequency is apparent (green shaded area). Without 
phase noise, the spectrum in blue would be detected, with an area corresponding to the 
actual phonon population. The dip from the phase noise background reduces the area 
under this curve, causing the area to be underestimated. In b the spectra from blue- 
side detection are presented. Here we note that the phase noise is amplified, causing 
a larger temperature. Though a and b are for large phonon occupation numbers, an 
asymmetry still appears due to classical laser phase noise. 

2.4-2. Effect on calibrated cooling beam thermometry Since direct photodetection is 
sensitive to the total field intensity and not the phase, laser light with phase noise 
detected without first passing through a system with frequency dependent transmission 
will fail to exhibit any fluctuations in excess of shot-noise. Ignoring the mechanical part 
of the system for the moment, a laser detuned from an optical cavity by A will, however, 
cause the phase noise of the laser near frequency A to appear as an increased noise floor 
level in the photodetected spectrum of the transmitted optical signal. The noise from the 
mechanical system undergoing random motion appears as a Loretzian peak on top of 



Laser noise in cavity- optomechanical cooling and thermometry 



13 



this noise floor in the photocurrent spectrum. Depending on the exact experimental 
geometry, i.e. whether the measurement is done on reflection or transmission and 
whether the probe laser is red- or blue-detuned from the cavity, the amplitude of this 
lorentzian signal above the raised noise floor, can be in excess or below that expected 
in an ideal measurement lacking laser phase noise. At very high relative noise levels, 
it is even possible for the signal peak to invert, and become a dip in the noise floor. 
Such an effect has been called noise squashing , and results from correlations in the 
input laser noise and the optically transduced mechanical motion noise due to radiation 
pressure fluctuations stemming from the same input laser noise. 

Figures |2^i and b display a model of the photocurrrent noise power spectrum for 
the transmitted laser field of an evanescently-coupled cavity-optomechanical geometry 
(Fig. [T^) for red (A = oj m ) and blue (A = — u m ) laser detuning from the optical cavity 
resonance, respectively. The shaded orange area denotes the part of the photocurrent 
noise power spectrum that is generated due to phase noise on the probe laser. It can be 
seen that the aforementioned interference effect, in the transmission geometry considered 
here, causes a dip (shaded green) in the phase noise background for red-side laser driving 
and a peak for blue-side laser driving} Note that, over a broader bandwidth than that 
shown in Fig. [2j the phase noise contribution to the photocurrent noise spectrum is 
modulated by the optical cavity lineshape, k » 7. 

Formally, the photocurrent noise power spectrum of Fig. [2] can be derived by 
considering the properties of a driven cavity-optomechancal system. The transmission 
and reflection of light by a laser driven optomechanical cavity has been of interest for 
a variety of switching and buffering applications [53], and displays physics analogous to 
electromagnetically induced transparency (EIT) [231 ESI [SSI S3 and electromagnetically 
induced absorption and amplification (EIA) [TTJ SHI E7] . Considering the evanescently 
coupled geometry of Fig. [T^i, the reflection coefficient of (weak) probe light at frequencies 
oj from an intense laser drive tone at frequency ul = u + A, is given by 

r±H = ^ ^ . (30) 

% A + u) + K/2 + -r. l %r /9 

Here reflection is into the backwards waveguide direction (a outr of Fig. [T^i). The 
transmission coefficient into the output channel in the forward waveguide direction is 
t±(co) = 1 + r ± (u). 

The photocurrent of the detected signal in the forward waveguide direction output 
due to laser phase noise present at the input is, 

1^(00) = t(u)a- m ^(u) + h.c.(-u), (31) 

which, taking into account the above relation between the transmission and reflection 
coefficients and the correlation between the positive and negative frequency components 
of the phase noise, yields in terms of the positive phase noise components only 

I (P (u)=r(u J )a^(u J )+h. C .(-u J ). (32) 

% In the case of detecting in the reflection channel for this geometry, the reverse is seen, with red-side 
laser detuning resulting in anti-squashing and blue-side detuning resulting in squashing. 
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Using eqns. (A. 3) and (26), the resulting photocurrent noise power spectrum in the 
output channel due to laser phase noise is then calculated to be: 

S v » = {\r{u)\ 2 6{u) + \ri-u)\ 2 6{-uj)) S EE {u). (33) 

Evaluating the reflection coefficient at laser detuning A = +u m yields for the full 
expression for a sideband-resolved system in the driven weak-coupling regime, 



T /^e\ 2 |7om| 7j = I^W 2 ! ^ 



k) 2 (w m + o;) 2 + (7/2) 2 ' 
The above expression in eqn. (34) is that plotted in Fig. [2j showing the laser phase 



noise contribution to the measured output noise power power spectrum. Adding this 



noise spectrum with detuning set to A = u m , to that found in eqn. (16) for the laser 
cooling beam output noise spectrum in the absence of classical laser noise, we find for 
the total transduced noise power spectral density near the mechanical resonance in the 
large cooperativity limit (ji « 7om), 

Sn(ou)\ A=Um = 1 + (—) + ^^—S bb (u; riinf), (35) 
\ K / 2k k 

where n m { = jiUb/j—Ken^^K and n b is mechanical mode occupancy in the absence of the 
cooling beam. This is the noise output power spectral density in experiments in which 
the laser cooling beam is also used for transduction/thermometry of the mechanical 
mode. Laser phase noise on the cooling beam, then, not only adds additional heating 



of the mechanical resonator (as captured by eqn. (29)), but the naively inferred phonon 



number represented by n- m { is also in error relative to the actual average mode occupancy, 

<n>jA=a> m ~ n ini = — n<p, (36) 

where we have assumed again the high cooperativity, deeply sideband resolved limit. As 
mentioned above, similar relations for heating and transduction error may be derived 
for other forms of classical noise, such as laser intensity noise and internal cavity noise. 

2.4-3. Effect on sideband asymmetry thermometry As indicated pictorially in Fig. [2]and 



mathematically in eqn. (34), the contribution to the detected photocurrent noise power 
spectrum due to laser phase noise takes on a different sign for red and blue detuned 
driving. In the evanescent coupling geometry considered, the former causes a dip in 
the phase noise background at the mechanical frequency resulting in noise squashing, 
while the latter leads to a reflection peak in the photodetected noise corresponding noise 



anti- squashing. For the sideband asymmetry thermometry described in Section 2.3, this 
classical noise asymmetry can mask the quantum asymmetry associated with zero-point 
fluctuations of the mechanical resonator, and thus must be carefully accounted for in 
such measurements. 
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From eqns. (16), (20), and (34) we find the inferred mechanical mode populations 



for probe measurements at detuning A = +uj. 



m j 



+ lin b 7om (k/2) 2 (K e \ 7; + 7om/2 

7+ 7+ (w m + uJ m ) 2 + (k/2) 2 V k / 7+ 
where n b is the average phonon occupancy of the mechanical resonator in the absence of 
the probe field. In the limit of low probe power (low cooperativity, C T = |7om|/7i « 1) 
applicable to measurements performed in Ref. [21] this simplifies to, 

«Mlc r «i = — — + — — n*. (38) 
7+ K 7+ 

Note that in this limit, averaging the two detuning measurements of nf ni results in 
a cancellation of the noise squashing for A = u m and the noise anti-squashing for 
A = — u m . This hints at a method of accurately determining the unperturbed phonon 
occupation number, i.e. n b = (j+n^ + j^n hii )/2j i , even in the presence of laser phase 
noise. 



For thermometry based upon the motional sideband asymmetry (Section 2.3), 
however, the effects of classical laser phase noise can not be so easily separated from the 
quantum noise asymmetry generated by the zero-point fluctuations of the mechanical 
resonator and the quantum back-action of the vacuum fluctuations of the probe laser. 
In the low cooperativity regime, the motional sideband asymmetry parameter described 
in Section 2.3 is modified to include the effects of probe laser phase noise, 

One way to sort out asymmetry effects related to classical laser phase noise, as is done 
below in Section 3.3, is to analyze the dependence of the measured asymmetry on the 
laser probe power. 



3. Experiment 

Below we describe, using the formalism presented above, several recent experiments [23l 
[21] involving the laser cooling and thermometry of GHz- frequency mechanical resonators 
formed within optomechanical crystals (OMCs) [13]. Additionally, we present new 
experimental results, on a more recent OMC cavity |58j system, being simultaneously 
cooled and probed using lasers with and without phase noise. As depicted in Fig. [3j 
OMCs are engineerable nanoscale structures that may be used to co-localize optical 
and acoustic (mechanical) waves. The particular device shown in Fig. [3J3 is a thin 
(220 nm) nanobeam formed from the Si device layer of a silicon-on-insulator wafer which 
supports a "breathing mode" mechanical resonance at a frequency of 5.1 GHz, and a 
high-Q optical resonance at an optical frequency of 194 THz (wavelength A » 1500 nm). 
Theoretical calculation of the optomechanical coupling strength between the breathing 
mode and co-localized optical mode yields a value of g /2iT = 860 kHz, due primarily in 
this design to the elasto-optic effect within Si. 
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Figure 3. a, Finitc-clcmcnt-method (FEM) numerical simulations of the electric field 
amplitude of the fundamental optical mode (theoretical frequency, lo /2tt = 194 THz) 
of the OMC cavity described in detail in Ref. [55]. b, FEM numerical simulation 
showing the displacement amplitude of the coupled breathing mechanical mode whose 
theoretical frequency is w m /27r = 5.7 GHz (measured at 5.1 GHz), c, Scanning electron 
micrograph of a fabricated version of the silicon nanobeam optomechanical crystal 
cavity. 





Table 1. 


Optomechanical 


crystal device and measurement parameters 




device 


A (nm) 


k/2-k 


K e /2ll 


g /2n 


U m /21T 


h/2tt 






(MHz) 


(MHz) 


(kHz) 


(GHz) 


(kHz) 


A [23] 


1537 


488 


65 


910 


3.68 


35 


B 


1547.3 


694 


97 


910 


5.1 


12.2 


C [21] (cooling) 1460 


390 


46 


960 


3.99 


43 


C [21] (read- 


1545 


1000 


300 


430 


3.99 


43 



out) 



Three different experiments are described below, each of which is performed with 
a slightly different OMC device and under slightly different optical coupling conditions. 
For reference below, in Table [T] we provide the important experimental parameters for all 
three devices, labelled A, B, and C. Device A was used in the experiment of Ref. [23] to 
cool a mechanical resonator close to its quantum mechanical ground-state, device B is a 
newly designed high frequency device [58J , and device C was studied in Ref. [24] as part 
of an experiment to measure the effects of quantum zero-point motion. As indicated, for 
device C, we have both a resonant optical mode used for cooling and a resonant optical 
mode used for read-out, with both modes coupled to the same mechanical mode. 

3.1. Experimental set-up 

The basic experimental set-up used for laser cooling and calibrated thermometry of 
the optomechanical crystal devices is shown in Fig. |4} A tunable external cavity 
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Figure 4. Laser cooling and calibrated thermometry experimental setup. 

As indicated by the separate boxed areas, the experimental setup consists of a tunable 
laser source used for cooling and mode thermometry, the OMC device under study 
which is placed in a continous-flow liquid helium cryostat, an optical pre-amplification 
stage used to amplify the cooling beam signal after transmission through the OMC 
cavity, a high-speed photodiode and electronic spectrum analyzer for measuring the 
photocurrent noise power spectral density of the transmitted cooling beam signal, 
and a modulation/lock-in scheme for probing the near-resonance optical reflection of 
the OMC cavity in the presence of the strong cooling beam tone (EIT-like reflection 
spectroscopy). The blue lines indicate the optical path for the cooling and mechanical 
mode thermometry measurement, while the dashed black lines indicate the alternative 
switched paths for calibrating the mode thermometry. The laser source is an external 
cavity semiconductor diode laser (New Focus Velocity laser, either model TLB-6328 or 
model TLB-6728), which can have its wavelength tuned over an approximately 60 nm 
wavelength span centered around A = 1550 nm. Other acronyms are: 2x2 optical switch 
(SW), variable optical attenunator (VOA), electro-optic modulator (EOM), lock- 
in detector (LI), erbium-doped fiber amplifier (EDFA), real-time spectrum analyzer 
(RSA), power meter (PM), optical detector (D), fiber polarization controller (FPC), 
optical wavemeter (A-meter), and RF signal generator (RF-SG). 



semiconductor laser source (New Focus Velocity series) is used both as the cooling laser 
beam, and upon transmission through the OMC cavity, to perform mode thermometry 
on the optically coupled mechanical breathing mode of the OMC device^ This 
experimental apparatus allows one to accurately set and measure the laser power at 
the input of the OMC cavity, and to calibrate the transduced mechanical noise power 
spectrum that is imprinted on the transmitted cooling laser beam. The series of micro- 
electro-mechanical 2x2 optical switches are used to switch the optical path repeatedly 
into a variety of different configurations, with very little variance (< 1%) in power 
levels. An optical fiber taper, formed by heating and stretching a single mode SMF28 
optical fiber down to a ~ 2 /zm diameter, is used to evanescently couple light into and 
out of the OMC device. This allows the use of fiber optics throughout the set-up for 
distribution of optical signals, providing a highly stable set-up in which measurments 

+ Measurements of the laser phase noise of each of the different laser models used in the measurements 



described below are presented in Appendix E 
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can be performed over days or even weeks. In order to provide a modicum of pre-cooling, 
the OMC devices are mounted into a continuous flow liquid helium cryostat, with an 
attainable temperature of ~ 6 K measured on the sample mount stagd^ A series of 
attocube piezo and slip-stick stages are used to position the fiber taper in the near-field 
of a chosen OMC device on the sample. Typically the fiber taper is aligned roughly 
parallel to the nanobeam OMC device, and placed in contact with the surrounding 
sample surface (but not the nanobeam itself) roughly 100 nm to one side of the OMC 
cavity. A Teflon Swagelok fitting with a pair of small diameter holes drilled in it are 
used to feed the optical fiber taper into and out of the crystat. 

During a laser cooling and thermometry run, a series of measurements are performed 
at each laser cooling power. These include: (i) the setting and stabilization of the 
cooling beam laser frequency to A = u m , (ii) calibration of the optical transmission, 
amplification, and detection, and (iii) measurement of the laser cooled mechanical noise 
power spectral density and noise background level. In the first step, the laser cooling 
beam frequency is set to a mechanical frequency red detuned of the optical cavity 
resonance (A = u m ) using near-resonance reflection spectroscopy of the driven cavity 
system. Such reflection spectroscopy, what we call EIT-like spectroscopy (see Ref. [17]). 
involves the use of a weak optical sideband of the intense laser cooling beam. The 
optical sideband is generated via the EOM, and is swept across a frequency span of 
A' = 1-8 GHz using the RF-SG. A small amplitude modulation is also applied to the 
optical sideband at a frequency of cjli = 100 kHz, allowing for lock-in (LI) detection 
of the reflected sideband signal from the OMC cavity. The cooling beam frequency 
is adjusted until the reflected sideband is aligned with the optical cavity resonance at 
a modulation frequency equal to the mechanical frequency, A' = w^JJJ The cooling 
beam laser frequency is then locked to within +5 MHz of this point using the optical 
wavemeter (A-meter). 

Calibration of the delivered laser cooling beam power to the input of the OMC 
cavity, and of the photodetected signal of the cooling beam transmission through the 
OMC cavity, are performed as described in detail in the Supplementary Information to 
Ref. [23J . In brief this involves measurment of the optical power at important (fiber taper 
input, fiber taper output, high-speed photoreceiver, etc.) points along the optical signal 
path shown in Fig. |4j Measurement of the transmitted power at detector (D3) is used as 
a reference, being constantly monitored using switch S W3 during the calibration process 

The temperature of the breathing mode is calibrated using the methods described in the main text 
and is used to extract the temperature TJ, of the mode's locally coupled bath. This bath temperature 
is typically higher, at Tj, m 10 — 20 K due to blackbody heating of the sample through the imaging 
viewport of the cryostat, and imperfect thermalization of the sample surface to the stages, 
jj Due to interference between light which is directly coupled into and out of the cavity, and light which is 
coupled into the cavity and interacts with the mechanical resonance, the corresponding reflected optical 
signal of the weak sideband has in addition to the normal broad Lorentzian feature of the optical cavity, 
a narrowband dip in the reflection. This reflection dip is analogous the transparency window in atomic 
EIT, and has a bandwidth of the laser cooled and damped mechanical resonance, with center frequency 
one mechanical frequency blue-shifted from the laser cooling beam frequency 
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and throughout subsequent measurements. Calibration also involves the measurement 
of the optically induced damping by the laser cooling beam, via measurement of the 
linewidth of the optically transduced thermal noise power spectrum of the breathing 
mode, for light sent in both directions through the fiber taper (enabled by the two 2x2 
optical switches at the input and output of the fiber taper, SW2 and SW3 of Fig. [1]). 
Asymmetries in the optically induced damping allow one to determine the optical loss 
before and after the OMC cavity in the fiber taper, which in combination with the 
total insertion loss of the fiber taper, provides an accurate estimate of the optical input 
power directly at the cavity. Calibration of the response of the high-speed photoreceiver 
(D2) and the optical amplifier (EDFA) is performed before every measurement point by 
applying an amplitude modulation on the input laser beam (using the EOM of Fig. [1]) 
of known amplitude at a frequency of the mechanical mode of the device under test, 
and recording the measured response in the photodetected (D2) noise spectrum. 

Finally, measurement of the laser-cooled mechanical noise spectrum requires three 
spectra for proper calibration: the mechanical noise spectrum taken with the cooling 
laser tuned to the optimal cooling point (A = u m ), a dark spectrum taken with the 
cooling laser far detuned from the cavity (A > Au m » k), and a background spectrum 
with the cooling laser blocked. Comparing the first and second spectra, we obtain the 
classical noise properties of the laser beam, while the relation between the first and 
second spectra provide us with the amplifier and detector shot noise levels inherent to 
the measurement apparatus. 



3.2. Laser cooling and cooling beam mode thermometry 

As described above, the mechanical mode occupancy can be determined by careful 
calibration of the optically-transduced mechanical noise power spectral density 
imprinted on the laser cooling beam. In the Supplementary Information of Ref. [23] 
we describe a method which can be used to both calibrate the laser cooled mechanical 
mode occupancy as well as to bound the effects due to classical laser noise. A plot of a 
typically measured noise spectrum in that experiment is shown in Fig. |5^l, taken under 
optimal cooling conditions with A = u m . The effect of classical laser noise on heating 
and thermometry can be discerned by the level of the noise background in the vicinity 
of the mechanical resonance peak. One can separate out additional classical laser noise 
from laser shot noise and detector noise in the measured noise background by making an 
additional noise measurement with the laser detuned by several mechanica frequencies 
from the optimal cooling detuning point, A > 4u> » k. This changes the measured 
classical laser noise in the vicinity of the mechanical resonance peak due to the filtering 



properties of the optical cavity. In the case of phase noise, as shown in eqn. (35), the 
background noise floor around the mechanical resonance drops by a factor proportional 
to 1 + (fte/ft) 2 ?^ when the laser detuning shifts from A = u m to A > 4u m . Similar 
relations can be derived in the case of laser intensity noise or intra-cavity noise (such as 
might result from thermo-refractive noise of the cavity). 
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Figure 5. a, Measured noise power spectral density near the mechanical resonance 
frequency of the transmitted optical cooling beam with detuning A = u> m . This data 
corresponds to the cooling point at m c = 10 intra-cavity photons of the cooling laser in 
Fig. 4b of Ref. [53], and appeared in the Supplementary Information of that work, b, 
Normalized background noise levels (see text) versus cooling beam intra-cavity photon 
number, again from the experiment reported in Ref. |23j . The green square point (□) 
was taken with the cooling laser beam filtered through a narrowband (50 MHz) optical 
filter prior to being input to the optomechanical cavity. 



Referring to the inset in Fig. |5^l, by taking the ratio between the difference in the 
two background noise measurements (ABG) and the noise peak of the phonon signal 
(A), we find for the case of laser phase noise 

ABG _ _ JH_ ^ 

A 2Kn inf 7oM 2Kn inf ' 
in the large cooperativity limit. Measurements from the ground-state cooling experiment 
of Ref. [23] are shown in Fig. [sJd for the normalized background level (ABG/A). These 
results show a normalized noise level less than 1.5% for all cooling beam powers (intra- 
cavity photon numbers). For the cavity coupling rates of the device measured in Ref. [23] 
(device A in Table [T]), this yields a bound on the laser phase noise photon number near 
mechanical resonance of n^,(a; m ) < 0.19. From eqn.(35) and Table [TJ for this level of 
laser noise the mechanical mode heating is [K e /2K)n^ < 0.012 quanta and the amount 
of noise squashing is {K e /K)n^ < 0.024 quanta, much less than the smallest inferred 
phonon occupancy n in f = 0.85 of the laser-cooled mode. As such, laser phase noise 
heating and squashing in these measurements were determined to be insignificant. A 
further verification of this is that no significant change in the noise background level or 
mechanical mode peak was detected with the cooling laser beam transmitted through an 
additional scanning Fabry-Perot filter (bandwidth 50 MHz), placed at the input to the 
optomechanical cavity. The normalized noise level for this filtered laser measurement is 
shown alongside the original data presented in Supplementary Information of Ref. [23] 
as a green square point in Fig. [5)d. 

Comparison of the above technique to more conventional laser phase noise 
measurements presented in Appendix E| can be made by converting the estimated noise 
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quanta into units of laser frequency noise. For the largest cooling beam intra-cavity 
photon number used in the experiment of Ref. [23] of n c % 2000 (corresponding to input 
laser power of \E \ 2 = 8.33 x 10 14 photons/s), the laser phase noise is bounded by (see 



eqn. (23)), S$ 



(n^/lEol 2 ) % 2.3 x 10~ 16 Hz" 1 at 3.68 GHz. For the Model 6328 laser 
used in this experiment, the corresponding laser frequency noise at 3.68 GHz is thus 
at most = ui^S (/,<(, & 1.2 x 10 5 rad 2 Hz, which is consistent with the measured laser 
frequency noise at this frequency of S^u * 7x 10 4 rad 2 Hz, shown in Fig. E2 1 using the 
calibrated Mach-Zender technique described in Appendix E This also explains why 
there is no obvious trend in the normalized noise level versus laser power of Fig. [5Jd; 
laser phase noise in these measurements is small and is masked by the uncertaintities 
in the measured background levels. 




ra/271 (GHz) co/2tu (GHz) 



Figure 6. a Broadband noise power spectral densities of the high frequency device B, 
measured using two different cooling beam lasers (detuning A = u m ), one posessing 
significant phase noise near the mechanical frequency (red curve; Model 6728), the 
other near quantum- limited (blue curve; Model 6328). b Zoom-in of the measured 
noise power spectral density around the cooled phonon peak at 5.1 GHz. 



In order to actually measure the deleterious effects of laser phase noise on the cooling 
and mechanical mode thermometry, we have also performed measurements of an OMC 
device with mechanical resonance at 5.1 GHz (device B in Table [TJ. From phase noise 
measurements over a wider span (1-7 GHz) of the Model 6328 and Model 6728 lasers 
studied in Appendix E, we found additional phase noise peaks at the second harmonic 
of the fundamental phase noise peaks shown in Fig. E2 For laser wavelengths near the 
optical resonance of device B (A = 1547.3 nm), this results in second harmonic phase 
noise peaks at 5 GHz and 6 GHz for the Model 6728 and Model 6328 lasers, respectively. 
The phase noise peak at (jj/2tt % 5 GHz of the Model 6728 laser was measured to have 
a corresponding frequency noise spectral density of S^ k = 7.5 x 10 6 rad 2 Hz, calibrated 
using the same techniques described in |Appendix E The measured phase noise of the 
Model 6328 laser is roughly 20 dB smaller than that of the Model 6728 laser at 5 GHz. 
The measurement was set up so that the cooling laser could be switched between the 
noisy Model 6728 laser and the quiet (in this frequency band) Model 6328 laser used 
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in the previously described ground-state cooling measurements above. Care was taken 
to match the laser power and polarization incident on the cavity for both lasers, thus 
enabling direct comparison of the measured results with differing levels of phase noise 
present. 

Figure [6] shows the results of such a measurement on device B, with red and blue 
curves corresponding to measurements with the noisy (Model 6728) and quiet (Model 
6328) lasers, respectively. Both of these measurements were taken for an optimal cooling 
detuning of A/2ir = u m /2ir = 5.1 GHz and an input laser power of \Eq\ 2 = 8-1 x 10 14 
photons/s (corresonding to an intracavity photon number of m c = 240). The effects 
of laser phase noise from the Model 6728 laser are immediately evident in the broad 
frequency scan of Fig. |6^, in which a broad noise peak can be seen around the mechanical 
resonance of interest at 5.1 GHz (there are in fact three narrow mechanical resonance 
lines visible in the data; the breathing mode mechanical resonance of interest which is 
most strongly coupled to the optical mode is the one at 5.1 GHz). Note that there is 
no such broad noise peak discernable for the blue curve of the quiet Model 6328 laser. 
A zoom-in of the measured mechanical resonance for both lasers is shown in Fig. [SJo. 

A fit to the area under the mechanical resonance yields an inferred phonon 
occupancy of n in f = 0.8 for the noisy Model 6728 laser and n m i = 1.2 for the quiet 
Model 6328 laser, can be estimated from the inferred phonon occupancy and the 
ratio of the mechanical resonance peak height (A) to the laser phase noise background 



(N BG ). From eqn. (35), = (N BG /A)(2K/K e )n- mi . For the red curve in Fig. |2p of the 
Model 6728 laser this yields % 3.2 noise quanta, corresponding to (K e /2/c)n^ 0.4 
phonons of heating and twice that of noise squashing, resulting in an actual mechanical 
mode occupancy of (n)± = 1.6 phonons (note that this is larger than the inferred 
phonon occupancy using the quiet Model 6328 laser by the 0.4 phase noise heating 
phonons, as one would expect). Putting in units of laser phase noise, we have 
from eqn. (g that S # = (n^/\E \ 2 ) * 4 x 10~ 15 Hz" 1 at 5.1 GHz for the Model 
6728 laser. From the broad spectrum in Fig. [2^i, the peak laser phase noise occurs at 
cu/27r = 5.05 GHz, and is approximately 1.65 times the phase noise at the mechanical 
resonance frequency, S , ^(cu/27r = 5.05 GHz) % 6.7 x 10~ 15 Hz -1 . This corresponds 
to a peak laser frequency noise of = u 2 S H % 6.8 x 10 6 rad 2 Hz, which is in good 
agreement with the independently calibrated peak laser frequency noise for this laser of 
7.5 x 10 6 rad 2 Hz. 



3.3. Sideband asymmetry mode thermometry 

In addition to the calibrated mode thermometry measurements described in the 
experiments above, we have also recently performed a form of self-calibrated 
measurement of the mechanical mode occupancy involving the measurement of 
the asymmetry in the motional sidebands generated by laser scattering from an 
optomechanical crystal resonator near its quantum ground-state [2U HI]. The device 
studied in this measurement is device C of Table [TJ and it has two optical modes coupled 
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to the same breathing mechanical mode of frequency u m /27T = 3.99 GHz. One of the 
optical cavity modes, resonant at wavelength A c = 1460 nm, is used to cool the breathing 
mechanical mode. The other optical mode, resonant at wavelength A r = 1545 nm, is 
used to read out the mechanical motion via the red (Stokes) and blue (anti-Stokes) 
scattered sidebands of a read-out laser beam near resonance with the read-out cavity 
mode. Further details of the measurement can be found in Ref. [24J, along with further 
discussion of the interpretation of the measurement results in Ref. |44j . Here we aim to 
re-present the data of Ref. |24j , in a form suitable for analysis of the effects of laser noise 



(primarily phase noise) as per the theoretical analysis described above in section 2.4.3 



We also compare the measured results with the expected noise effects from the calibrated 



laser phase noise of the read-out laser presented in |Appendix E 



Figure shows a plot of the inferred phonon occupancy of the breathing mechanical 
mode at 3.99 GHz versus the mechanical linewidth of the beathing mechanical mode 
(7 C ) for each cooling laser beam power. The cooling laser, a Model 6326 Velocity laser 
in this case, is tuned to a mechanical frequency red of the fundamental optical mode at 
A c = 1460 nm (A c = u m ). The inferred phonon occupancy in the absence of a read-out 
beam, n bc is measured using a read-out laser (the same Model 6328 Velocity laser used 
in the ground-state cooling experiments of Ref. [23]) that is tuned near resonance of the 
second-order optical mode of the OMC cavity at A r = 1545 nm. This occupancy {n bc ) 
can be thought of as the occupation number of a new effective thermal bath coupled 
to the mechanics consisting of a combination of the intrinsic mechanical damping and 
optical damping from the cooling beam, nbc is measured by taking the average of the 
calibrated phonon occupancies measured for read-out beam detunings of A r = ±u m ; 
n bc = l/2(7 + n + + 7_n_)/7, where n + and n_ are the inferred phonon occupancies from 
the measured read-out beam photocurrent spectrum and 7+ and 7_ are the measured 
mechanical linewidths for A r = u m and A, r = — u m , respectively. 7 is the mechanical 
damping rate of the total system in the absence of the a read-out beam and is calculated 
from the measured probe beam spectra using 7 = 1/2(7+ +7-)- At each of the cooling 
beam powers, the read-out laser beam power is adjusted such that the cooperativity of 
the read-out beam is substantially below unity (see Ref. [24J ) . As such, using this two- 
mode approach, one high power for cooling and one low power for read-out, separates 
the noise heating and noise squashing effects, and substantially reduces any effects due 
to laser phase noise on the mode thermometry. 

In Fig. [7)3 we plot the measured asymmetry in the motional sidebands, 77 in 
eqn. (21 ), of the read-out laser beam versus the inferred laser-cooled mode occupancy of 
Fig. [7^,. The red solid curve of Fig. [7]d also plots the expected sideband asymmetry from 
the linear fit of the laser-cooled mode occupancy versus mechanical damping (dashed 
black-curve of Fig. [7k). The asymmetry is quite pronounced at lower mode occupancies 



as expected, and matches well the theoretical sideband asymmetry curve of eqn. (21) 
due to the behavior of the zero-point fluctations of the mechanical mode and interference 
with the quantum back-action of the read-out laser shot-noise [2U 133] . As discussed in 



Section 2.4.3 such sideband asymmetry can also arise from classical laser noise on the 
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Figure 7. a and b show laser cooling and motional sideband asymmetry data from 
the experiment reported in Ref. [24 . In this experiment a separate cooling beam and 
probe beam were used, each coupled to the localized breathing mechanical mode of a 
Si nanobeam resonator at frequency uj m /2n = 3.99 GHz (device C of Table [l]). In a we 
plot the measured optically enhanced mechanical damping rate (7) versus the inferred 
laser cooled mode occupancy of the breathing mechanical resonator (ri& c ). Here nt> c 
represents the breathing mode occupancy due to laser cooling in the absence of the 
probe beam. It is computed as rl& c = + Ti^f)/27, from the inferred mode 

occupancies (n.J f ) and mechanical lincwidths (7+ and 7 = (7+ + 7_)/2) measured 
with the weak probe beam at detunings A r = +uj m . As plotted in b, the weak probe 
beam signals for detunings A = ±w m are also used to calculate 77 of eqn.(21 ). The red 
line in b represents the expected theoretical asymmetry (77 = l/n{, c ) due to zero-point 
motion of the mechanical system in the absense of phase-noise effects. In c we also 
plot the measured asymmetry versus the probe read-out power as represented by the 
read-out intra-cavity photon number m r . We have highlighted two subsets of data 
within the full set of data, represented by green and red circles. For the green circle 
data points the measured asymmetry is roughly constant at rj = 0.23 + 0.02 while the 
probe beam power is varied over a factor of nearly two. For the red circle data, the 
probe beam power is held fixed at n r = 40 and the asymmetry is measured to range 
over a factor of three. 



red-out laser beam (see eqn.(39)). Owing to the fact that the read-out laser beam is 
of substantially lower power than that of the cooling laser beam in these experiments, 
the amount of classical laser phase noise is expected to be small. Also, the read-out 
laser beam power was only weakly correlated with the cooling beam power (and thus 
laser-cooled mode occupancy), making it highly unlikely that the measured motional 
sideband asymmetry curve of Fig. [7]d is due to classical laser noise on the read-out beam. 
Verification of this is shown in Fig. [7]3, in which a scatter plot of the measured sideband 
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asymmetry is plotted versus the read-out beam laser power (intra-cavity photon number, 
n r ). We have highlighted in this set of data two groups of points. The red circles 
correspond to a range of data points in which the read-out beam power is varied over 
a factor of nearly two (n r % 30-60), whereas the measured asymmetry varies by less 
than +10% of its nominal value (r] = 0.23 + 0.02). The green circle points show data 
in which the read-out beam power was held fixed (n r » 40) as the cooling beam power 
was varied. For these data points the measured asymmetry is seen to vary by a factor 
of almost three (rj % 0.15-0.45). Both of these groups of data points serve to indicate a 
lack of correlation between the measured motional sideband asymmetry and the read- 
out beam power, strongly ruling out laser phsae noise as the source of the measured 
asymmetry in the experiments of Ref . [24J . 

Finally, one can estimate the magnitude of the effects of laser phase noise of 
the read-out laser on the measured sideband asymmetry in these measurements by 
comparing to the calibrated laser phase noise presented in |Appendix E From Fig. |E2^ i, 
the laser frequency noise of the Model 6328 laser used in the read-out of the sideband 
asymmetry experiment has a value of S^u, % 5 x 10 4 rad 2 Hz at a frequency of 
u)/27r = 3.99 GHz (the breathing mechanical mode frequency of device C) and for a laser 
wavelength near A r = 1545 nm. From Fig. [7}^, the read-out beam laser power is at most 
\Eo\r = n r /(ft e /2u^) < 4 x 10 13 photons/s. The maximum read-out laser phase noise 
quanta in these measurements is then bounded by = (S^/m^IEoH ^ 0.003 quanta. 
The corresponding laser noise heating of the mechanical resonator by the weak read- 
out beam is (C r /1 + C r )(/c e /2«;)n^, which for the read-out beam cooperativity of these 
experiments (C r < 0.1) is at most a negligible 4.3 x 10~ 5 phonons. Laser phase noise 
squashing and anti-squashing results in the modified motional sideband asymmetry given 
in eqn. (29 ) for a low-cooperativity read-out beam. The correction factor to the quantum 
asymmetry due to classical laser phase noise is given by (l + 2K e n^//c)/(l — K e n^/Kn^) — 1, 
which for the measurement of Ref. [21] is smaller than 0.2%. We should also note that 
from the measured phase noise of the Model 6326 laser used to cool the mechanical 
mode in these experiments (see Fig. E2:,f), the estimated laser phase noise heating 
from the cooling beam is less than n$ % 0.04 phonons at the largest cooling beam 
powers (n c = 330 intra-cavity photons), indicating that the minimum measured phonon 
occupancy of 2.6 is also not limited by laser noise in this experiment. 



4. Conclusions 



We have presented a unified analysis of the effects of optical noise, due to both laser 
phase noise and quantum fluctuations of the electromagnetic field, on thermometry and 
cooling in recent optomechanical experiments. By doing so, we provide a systematic 
means by which the presence of classical laser noise can be detected, and taken 
into account, and we further rule out any role which it may play in our recent 
experiments [23| 124"] . This result is particularly useful for experimenters working on 
gigahertz nanomechanical resonators, since it allows, with some care, for investigations 
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to probe qauntum optomechanics in the telecom wavelength bands with cheaper, and 
more readily commercially available diode lasers. 
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Appendix A. Definitions 

Fourier Transforms are defined for operators and variables in the symmetric manner 

1 f 00 

A(t) = -j= du e- lut A(u), 



00 

oo 

iuit 



A(u) = -j= dt e lUJt A{t). (A.l) 
Spectral densities are defined as 

Saa(co)=\ dr e™ T (ti{t + r)A(t)). (A.2) 

J— oo 

and symmetrized as Saa (w) = ^(Saa(w) + Saa(— Here the angular brackets denote 
expectation values as defined in quantum mechanics {A) = Tr(Ap). When classical 
stochastic processes are placed in angular brackets, we refer to a classical ensemble 
average. In Fourier domain, the various system operators are written in terms of the 
bath noise operators, and therefore knowledge of the expectation values of form (A(lu)) 
and (B(cu)A(uj')) where A and B are bath field operators is sufficient to calculate spectral 
densities. These correlations are known from the density matrix of the baths, which in 
this paper are assumed to be in either a vacuum or thermal state. The Hermitian 
conjugate of operator A{t) is A^(t), and has a Fourier transform denoted as 

it (a;) = -L f dt e**A\t), 

V27T J-oo 

leading to (A(u>)\ = A^(-uj). The spectral density may be written also as 

S AA (u) = f du/ <it(w)i(a/)>. (A.3) 

J— oo 
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Appendix B. Mechanical resonator spectral density 

We use mainly the input-output formalism to derive the various spectra which are being 
measured. For the case of the mechanical resonator, the quantum Langevin equation is 
given by 

S(t) = - (ioj m + bit) ~ V%L(t), (B.l) 



2 

With correlation functions for the input noise 

(bl a (t)L(t')) = n b 6{t-t'), (B.2) 
(b- m (t)bi(t'))=(n b + l)6(t-t'), (B.3) 

where n b is the occupancy of the thermal bath connected to the mechanical resonator. 
Inherent in the correlation functions above is the assumption that the mechanical 
resonance bandwidth is very small compared to its resonance frequency (high mechanical 
Q-factor), such that the bath occupation can be taken as a single number n b {u m ). 
Additionally, we find for the Fourier transform of the input noise operators, 

(b\ n (u)b in (uj')) = n b 6(u + u'), (B.4) 
(b- m (u)bl(u')) = (n b + l)5(u + <J). (B.5) 



By solving the Fourier transform of equation (B.l), we find that the mechanical mode 



annihilation operator will be deterimined by the input noise as 

S( W ) = , -^ L{UJ) , . (B.6) 
This can easily be used to calculate the spectral density, and we find, 

rco 

S bb (u) = dre WT (b\r)b) 



J — oo 

rco i poo rco 

dr e luJT — du" du' 

J -co ^ J -co J -co 

x ({b{u")yh{u'))e tLU " T 



( a;m + a; )2 + ( 7 . / 2)2- ^ B - 7 ) 

We will sometimes denote this function as S bb (u; n). This spectra density can be thought 
to represent the ability of the mechanical system to emit energy. A similar expression 
can be found for the creation operators: 

£W(w) = t 7 '^ 2 , 7 /o xn (B.8) 

{u m - u) 2 + (7i/2r 

This spectra density can be thought to represent the ability of the mechanical system 
to absorb energy. Finally, for the position operator, x = x zp f(b + ¥), we find [6] 

S xx (u) = x\ f (S bb (u) + S b tm(u)) 



„2 



X zp { I ~ ; 7CZT7C + 



ji-n 7i(n + 1) 



(u m + u) 2 + ( 7i /2) 2 (u m - u) 2 + ( 7i /2) 2 
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Appendix C. Quantum noise squashing 

Assuming perfect detection (n e /2 = k) we find that the normalized heterodyne current 
is given by 

2C 

I(t) = -ia- m (t) + idi(t) + —(b(t) + &t(t)) (C.l) 



Taking the autocorrelation of the detected current, 
</(r)/> = (a in (r)al) + \' ' 

K X zpf 



2,0 \aUrp) - <6(r)aL>) (C.2) 



K 

we find Sii(oj) by taking the Fourier transform of the above expression. The first term, 
due to the fact that the noise is delta correlated gives a constant noise floor. The second 
term can be thought of as a measurement of position, and we see that the rate at which 
information is gathered about the system is 4|G| 2 /k, i.e. the optomechanical damping 
rate, and back-action. The cross-correlation terms are calculated as such: 

dr e lWT (a- m {rp) = 



I 

1 

2^ 



00 

oo r ->- r 

dr e iulT 

-co J— oo J — co 



rCO rCO 

I du' I du" e- iuJ ' T (d in (u')(b(u")y) (C.3) 

J — CO J — CO 

Using the back-action modified mechanical fluctuation operator shown in equation ([5]), 

b(U) = — r ; 1 r — (C.4) 

i\u m — u) + 7/2 y/K i(u m — u) + 7/2 



and the properties of vacuum fluctuation operators, we find 

dr e^(a ia (rp) = ^— 1 (C.5) 

-oo V K "H u m 7/2 

00 InC 1 

dr e^(b(r)al) = + — , -. (C.6) 

H i{u m - u) + j/2 



I 



From here we calculate the spectral density of the heterodyne signal, 

S n (u) = 1 + 

4|G| 2 / 7n 7(n + 1) 



+ — -, ^ — t-^tk + 



k V k + w) 2 + (7/2) 2 (u m - u) 2 + (7/2) 2 

4|G| 2 7 
« (w m - w) 2 + (7/2) 



2 



= 1 + ^S bb {u). 

K 

Appendix D. Scattering matrix elements 

By algebraic manipulation of the Heisenberg equations of motion in Fourier domain, 
and an additional input-output boundary condition a out = a in + ^n e /2d, we arrive at 
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scattering relations 

a utM|A=-w m ^ t(u; A)a in (u)) + n opt (oj; A)a inji (a;) + s 12 (uj; A)b\ n (uj)(D.l) 

and 

o utMlA=u> m » A)a in (a;) + n opt (a;; A)a in)i (u;) + Si 2 (a;; A)6 in (a;), (D.2) 

for blue- and red-side laser pumping, respectively. In the driven weak coupling regime 
(7om « k), these scattering coefficients have simple algebaric forms which are presented 
below. 



Appendix D.l. Red-side driving: A = u m 

The red-side scattering matrix elements for values of u about the mechanical frequency 
(00 — u m « k) are: 

±1 a \ 1 K e , 70M «e 1 m q \ 

t(w;A = w m ) = l + — 77 r- D.3 

< A ^ /2^7/ |7om|/2 \ 

"opt w; A = w m ) = \ — — - 1 , D.4 

V /c \i(<4 - w) + 7/2 / 




si 2 (w; A = w m ) = * / — — (D.5) 



and 

Smfll) 1 A = f;).A = * / 

2Ki(u! m — id) + 7/2 

Appendix D.2. Blue-side driving: A = — w m 

The blue-side scattering matrix elements for values of u about the mechanical frequency 
(u + uj m « k) are: 

t(u; A = -u m ) = 1 - ^ - Jto^ 1 (D . 6) 

v ' ; « 2k -z(w m + w) + 7/2 v ; 



2 « / K e / |7om|/2 
k 2 \— z(w m + w) + 7/2 



n opt (a;; A = -u m ) = -\ — — ( — — — ; -3 ; - /o + 1 ) (D.7) 



/ a \ K e IV7< 70M m s x 

312(0;; A = -U m ) = \ ~ Tf ; r- (D.8) 

V 2k — «(u; m + lo) + 7/2 



Appendix E. Laser phase noise measurement 



The constraints of silicon nanofabricated devices, i.e. their single-mode nature, and 
their spread in parameters such as optical cavity frequency caused by fabrication 
imperfections, mean that wideband tunable external cavity diode lasers (ECDL) are 
invaluable for experiments in optomechanics. As such, they have been used extensively 
by our group as well as others [59j EH ETJ , for a variety of systems operating at different 
mechanical frequencies. Phase and frequency noise have been of concern in these types 
of lasers, and it is therefore important to evaluate the laser noise properties. 
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Here we present a measurement of See{u) based on a measurement of the phase 
noise spectral density S^^oS) and relation (23). A Mach-Zehnder intereferomoter 
(MZI) is an optical component where the detected intensity of the transmitted light 
is dependent on the frequency of the light, and given by T{u) = 1 + sin(w/co>FSR), where 
c^fsr is the free spectral range of the interferometer. For the moment, we linearize 
this relation as, T{uj) = 1 + cj/wfsr (the full relation is used for the presented data). 
Assuming that the laser frequency ojl = is a stochastic process representing 

the instantaneous frequency of the laser light, the detected intensity will be given by 
I(t) = N + NcuL(t)/u)FSR + ^sn(^)- The term N is the average flux of photons incident 
on the detector, while the last term nsnit), is the white shot-noise of the laser with 
amplitude proportional to A 1//2 . The spectral density of the detected signal is then 
given by, 

Sn(co) = N + —j—Suw(w)- (E.l) 

^FSR 

Therefore, using the shot noise as a reference level, we can define a signal-to-noise ratio 

as, 

SNR=^-S uu {u), (E.2) 

^FSR 

where the "signal" is the classical laser frequency noise and the "noise" is shot-noise. 
This shot-noise based calibration method can be performed by a set-up such as the one 
shown in Figure |Elb . Such measurement of shot-noise however depends on knowledge 
of the detector quantum efficiency 

An alternate calibration method, free of any need to measure absolute powers and 
detector efficiencies, or to model the optical component transducing the phase noise 
(here the MZI), is shown schematically in Figure Eli. A phase modulator is used to 
generate a tone with a large modulation index /3 at the frequency of interest. This 
unitless modulation index is directly obtained by measuring the power of the generated 
sidebands, using a scanning Fabry-Perot filter in our case. These powers are related to 
each other by the appropriate Bessel function with argument (3. As described below 
in Appendix E.l, measuring the tone passing through the MZI provides us with a 



calibration of S uu (u). 

Both methods give nominally the same results for the calibration of laser frequency 
noise, which we plot in Fig. E2i-c for three different New Focus Velocity series external- 
cavity semiconductor diode lasers used in the work presented here and in recent cavity- 
optomechanical experiments with GHz-frequency OMCs [2"3"l 12%] . In these plots, the 
detector noise at zero input power (typically denoted as noise-equivalent power (NEP)) 
and the laser shot noise, are subtracted from raw measured noise spectra. Spectra 



showing the various stages of the processing are shown in Fig. [E3 
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Figure El. a The experimental setup used for measurement and calibration of 
laser phase noise. Light from the laser under test is sent into an electro-optical phase 
modulator (0EOM), generating a pure phase modulation which is used to calibrate the 
system response. Most of the light is then taken from the 90% end of a beam splitter, 
and sent through an all-fiber imbalanced (AL ~ 2 cm) Mach-Zehnder interferometer 
(MZI), with a free spectral range of 20.1 GHz (measured using calibrated Toptica 
wavemeter). The MZI converts frequency fluctuations into intensity fluctuations, which 
are detected on a high speed photodetector (New Focus 1544-B). The DC output 
of the photodetector is used to perform a low-frequency lock of the laser frequency 
to the mid-point of the MZI sinusoidal transfer function. The 10% split-off signal 
is sent through a scannable optical filter (bandwidth 50MHz), which allows one to 
independently measure the modulation index of the phase modulated calibration tone 
by sweeping the filter across the carrier and sidebands, comparing their optical power. 
Fiber polarization controllers (FPCs) are used to adjust the optical field polarization, 
b For an alternate calibration of the laser phase noise, the shot noise level of the 
detected laser signal is measured. This is accomplished by sending the laser through 
a narrow optical filter (50 MHz bandwidth), followed by a variable optical attenuator 
(VOA), to obtain the shot noise level of the photodetected signal at GHz frequencies. 
The optical filter is used to remove classical laser noise at these frequencies. 



Appendix E.l. Phase noise calibration 



Here we describe the phase noise calibration method (with the setup shown in Fig. El i) 
used to characterize the lasers. An electro-optic phase modulator is used to generate 
sidebands on the optical laser signal at ul, by creating a phase modulation of 

<f>(t) = (3cos(u c t). (E.3) 

The ratio of the power between the carrier at ujl and the first order sideband is given 
by 



Pi 
P 



MP) 



(E.4) 



MP) 

Using a scanning Fabry-Perot filter with a bandwidth much smaller than u c , we select 
out each sideband individually, and measure the powers Po.i i n the carrier and sidebands 
to obtain a value for (5. 

The frequency noise spectrum for a known modulation <fi(t) can be calculated 
from the Fourier transform of the autocorrelation function (0(t)0(t + r)). For the 
case of sinusiodal phase modulation we have (<fi(t)<fi(t + r)) = u 2 (3 2 cos(o; c r)/2, with 
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Figure E2. a, b and c are the calibrated frequency noise for Lasers 1 (Model 
6328, Serial #286), 2 (Model 6728), and 3 (Model 6326, Serial #19), respectively. In a, 
b the different color curves are for measurements made at different laser wavelenghts, 
with {blue, green, red, cyan} corresponding to A ={1520, 1537, 1550, 1570}nm, 
respectively. Laser 3 is a 1400nm band laser which was only operated at A = 1460 nm. 
The plots shown in d, e show the corresponding level of laser phase noise heating 
{firj,) in units of phonons for the device parameters and maximum cooling laser power 
(n c = 2000 intra-cavity cooling beam photons, K e /2ir = 65 MHz, k/2tt = 488 MHz, 
ui m /2ir = 3.68 GHz) of Ref. [23]- In f we show the same for the device parameters 
and maximum cooling beam power (n c = 330, K e /2w = 45 MHz, h/2-k = 300 MHz, 
w m /27r = 3.99 GHz) of Ref. [22]. 



a corresponding frequency noise power spectral density of (in units of rad Hz), 

S c fM = ^ - We ) + 5{u + Wc )) . (E.5) 

By comparing the raw measured noise of the laser frequency noise (S meas (co)) to 
that of the raw measured noise in the calibration signal peak (S^^lj)), one can 
calibrate the measured laser frequency noise in units of rad 2 Hz. Specifically, we have 
that A(u e ) Z-Au S%L(w)du = ttw 2 /3 2 /2, where A(oo c ) is the conversion coefficient 
(at frequency lu c ) between measured electrical noise power denity and (symmetrized) 
frequency noise power density for our experimental apparatus. The corresponding laser 
frequency noise at u c can then be related to the measured electrical noise and the noise 
power in the phase modulation calibration tone as, 

SuuM = ^ J mUUc) ■ (E.6) 



In order to generate the calibrated laser frequency noise spectra of Fig. E2 we measured 
the conversion coefficient, A(u c ), at ~ 50 MHz intervals across the entire frequency span 
of the measurement. 




Figure E3. a, Measured total noise spectrum (red curve), along with the shot noise 
(blue curve) and detector NEP (yellow curve), used to generate the laser frequency 
noise spectrum of Fig. E2 1 (only 1520 nm wavelength shown). The detector NEP is 
measured by simply blocking the laser (zero optical power) and measuring the noise 
spectrum (note that the detector NEP is well above the intrinsic noise level of the 
spectrum analyzer). The shot-noise is measured independently for an identical laser 
power and wavelength using the set-up of Fig. |Elb . The green line corresponding 
to the sum of average NEP and shot noise levels. The black curve is the measured 
signal with the shot noise and NEP background levels subtracted, and represents the 
component of the measured noise which we attribute to laser frequency noise in Fig. 



E2i. b, Measured total noise spectrum (red curve) and detector NEP (yellow curve) 



used to generate the laser frequency noise spectrum of Fig. [E2b. The shot-noise level 
(not shown) for this measurement is far below the NEP level. The black curve is the 
measured signal with the NEP subtracted, and is the resulting laser frequency noise 
spectrum plotted in Fig. |E2b. 
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